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We report self-consistent ab initio calculations of structural 
and electronic properties for the B- and Mg-terminated MgB2 
(0001) surfaces. We employ ultra-soft pseudopotentials and 
plane wave basis sets within the generalized gradient approx- 
imation. The surface relaxations are found to be small for 
both B- and Mg-terminated surfaces. For the B-terminated 
surface, both B a and tt surface bands appear, while only 
one B TT surface band exists near the Fermi level for the Mg- 
terminated surface. The superconductivity of the MgB2 sur- 
faces is discussed. The work function is predicted to be 5.95 
and 4.25 eV for the B- and Mg-terminated surfaces respec- 
tively. The simulated scanning tunneling microscopy images 
of the surfaces are not sensitive to the sign and value of the 
bias voltages, but depend strongly on the tip-sample distance. 
An image reversal is predicted for the Mg-terminated surface. 



The surprising discovery of superconductivity at 
39K in MgB2 has attracted great scientific interest. 
The isotope effect experiments [^[^ indicate that MgB2 
may be a BCS phonon-mediated superconductor, with 
Tc above the commonly accepted limits for phonon- 
assisted superconductivity. Band structure calculations 
show that the Mg donates substantially its 3s elec- 
trons to the B layer, and that the Fermi surface is derived 
mainly from B orbitals with two-dimensional character. 
The superconductivity is believed to be due to the strong 
coupling of the holes of the B a bands to the in-plane B 
phonon lattice vibrations (E2g modes) 

Being a simple, low-cost and high-performance binary 
intermetallic compound, MgB2 is a very promising can- 
didate for superconducting device as well as large-scale 
applications. Experimentally, scanning tunneling spec- 
troscopy measurements of the surface of superconducting 
MgB2 have been reported and the high-quality c- 

axis-oriented epitaxial MgB2 thin films have been grown 
successfully on AI2O3 substrates by using a pulsed laser 
deposition technique However, there has not been a 
theoretical study for the MgB2 surfaces so far. In this 
Letter, we report a first-principles study of MgB2 (0001) 
surfaces. Both B- and Mg-terminated surfaces are in- 
vestigated. We show that while the relaxations of the 
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surface layers are small for both surfaces, their electronic 
structures are quite different. 

Our calculations were carried out within the general- 
ized gradient approximation (GGA) j|] using ultrasoft 
pseudopotentials and plane wave basis sets A 
1x1x15 supercell geometry was used to model the sur- 
faces in which there are 15 atomic layers and 15 layers of 
vacuum (hereafter, 15-layer slab model) . The plane wave 
cutoff is 257.2 eV for structural optimizations and 321.5 
eV for static electronic structure calculations. Brillouin 
zone integrations were performed on a grid of 13 x 13 x 1 
Monkhorst-Pack special points. During the struc- 
tural optimizations, we fixed the central 3 atomic layers 
in the bulk configuration and allowed all other atoms in 
the supercell to move until all forces vanish within 0.02 
eV/A. We have checked some our calculations with a 17- 
layer slab model, and found the 15-layer slab model is 
sufhcient for our studies. 

The bulk MgB2 consists of alternating hexagonal lay- 
ers of Mg atoms and boron honeycomb layers. Our bulk 
calculations give the lattice constants of a = 3.067A, 
c — 3.515A, which are in good agreement with experi- 
mental and other theoretical results [Hl3|. We used this 
bulk geometry as the initial geometry for the surface op- 
timizations. Table 1 lists the results for the surface relax- 
ations, given in terms of the change of the interlayer 
distance in percent of the distance in the bulk. From 
this table, one can see the relaxations for both the B- 
and Mg-terminated surfaces are rather small (less than 
4%). The relaxation affects up to five atomic layers for 
the B-terminated surface, while it mainly localizes on the 
first two layers for the Mg-terminated surface. For the B- 
and Mg-terminated surfaces, an inward relaxation of the 
top layer by -2.1% and -3.7% and an outward relaxation 
of the second layer by 2.0% and 1.2% are observed. The 
energy gains due to the relaxation are 2.31 and 2.38 eV 
for the B- and Mg-terminated surfaces respectively. 

The surface electronic structure of the B- and Mg- 
terminated MgB2 (0001) surfaces at the optimized ge- 
ometries is shown in Fig.l together with the projected 
band structure (PBS) of the bulk MgB2. From Fig. 1, 
we can see that for the B-terminated surface, there is one 
surface band along each of the bulk bands. The disper- 
sion of these surface bands is very similar to that of the 
bulk bands, and they can be sorted into the bonding a 
(spa;Pj,) and nonbonding tt (pz) bands of the surface B 
atoms accordingly. For the Mg-terminated surface, how- 
ever, there is only one surface band outside the PBS. This 
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band mainly results from the orbitals of B atoms in 
the second surface layer. 

The above results can be understood by plotting the 
charge density distribution. The charge density contour 
maps for the x-z plane and x-y plane at B layers are 
shown in Fig. 2. In good agreement with other calcula- 
tions for the bulk, we found that electrons are transferred 
form Mg orbitals to B orbitals, but they are not well lo- 
calized at the B sites, instead they are distributed over 
the whole crystal. The MgB2 is essentially metallic boron 
held together by covalent B-B and ionic B-Mg bonding. 
At the B-terminated surface, both in-plane (spa;Pj,) and 
out-plane (p^ ) electrons of the surface B layer are slightly 
redistributed due to lack of one nearest neighbor (NN) 
Mg layer. This kind of electronic redistribution results 
in the appearance of the B a and tt surface bands. At 
the Mg-terminated surface, however, the second surface 
B layer has two NN Mg layers as in bulk, and only its 
out- plane (p^) electrons are perturbed by the electron 
redistribution of the surface Mg layer, leading to a B tt 
surface band. 

Figure 3 shows the total densities of states (DOS) 
and layer resolved projected DOS for the B- and Mg- 
terminated MgB2 surfaces together with those for the 
bulk. While the total DOS curves are quite similar over- 
all, there are differences in the layer resolved project 
DOS. The surface-specific features in DOS disappear 
slowly for the B-terminated surface, while they vanish 
rapidly for the Mg-terminated surface. For the Mg- 
terminated surface, those features are mainly localized 
in the (Mg, B) surface layers, the DOS of the first (Mg, 
B) sub-surface layers are already virtually identical with 
those of the bulk. Even in the second (B, Mg) sub-surface 
layers for the B-terminated surface, however, we still can 
see those features. The total DOS at the Fermi level 
(iV(0)) is calculated to be 0.42 and 0.37 states/(eVx 
B atom) for the B- and Mg-terminated surfaces respec- 
tively. Both are larger than the A^(0) of the bulk (0.36 
states/ (eVx B atom)). 

To discuss qualitatively the superconductivity of the 
MgB2 surfaces, we calculated the layer resolved zone- 
center (in-plane displacements of borons) phonon 
frequencies. They are 596 cm"-'^ and 479 cm^^ for the 
first surface B layer of the B- and Mg-terminated sur- 
faces respectively. These values are bigger than that of 
the bulk-like B layer ( 476 cm~^). From our A^(0) and 
E2g results, we expect that the superconductivity of the 
MgB2 surfaces is somewhat strengthened. This conclu- 
sion is contrary to the conjecture of a weakened surface 
layer The conjecture was resulted from that the 

measured superconducting energy gap is less than the 
weak-coupling BCS value of 5.9 meV for the bulk Tc of 
39 K 1^,^,0. However, we notice that the measurement 
might be done on chemically modified surfaces, not on 
the ideal surfaces. 

We also calculated the work function of the MgB2 sur- 



faces in terms of the difference between the self consis- 
tent potential in the vacuum and the fermi level. As 
listed in Table 1, the B- termnated surface gives a work 
function of 5.95 eV, while the Mg- termnated surface is 
4.25 eV. This result indicates that the electrons on the B- 
terminated surface are more stable than Mg- terminated 
surface. 

The scanning tunneling microscopy (STM) is a widely 
used tool in surface research. Here we simulated the 
STM images for the MgB2 surfaces using the Tersoff 
and Hamann theory [ p^ . By assuming an asymptoti- 
cally spherical tip and taking the limit of small applied 
bias voltage, the tunneling current was reduced to the 
local density of states (LDOS) of the surface at the point 
probed by the tip: 

I{x, y) cc p{x, y) = Y. - (1) 

i 

We found while the simulated images are not sensitive to 
the sign and value of the bais voltages, they depend on 
strongly the distance between the tip and sample. Fig. 
4 shows some simulated images with different tip-sample 
distances at -1-0.5 V. From these figures, one can see that 
at the Mg-terminated surface, the STM images show a 
"bright-dark" transition. At a short distance, the bright 
area of the image is corresponding to Mg sites in the sur- 
face layer. At a long distance, the bright area is changed 
to the B sites of the second surface layer. This rever- 
sal contrast of the STM images at different tip heights 
can be understood by Eq.(l). In Eq.(l), the tunneling 
current is determined by the space and energy windows 
simultaneously. At the Mg-terminated surface, the en- 
ergy contributions to the tunneling current mainly come 
from the B layer while the surface Mg layer gives the most 
space contributions when the tip is close to the surface. 

In summary, we have studied the structural and elec- 
tronic properties of the B- and Mg-terminated MgB2 
(0001) surfaces. The surface relaxations are found to be 
small. The surface bands are all along the bulk bands. 
For B-terminated surface, both po- and Ptt surface bands 
are exist, but for Mg- terminated there is only a p^r sur- 
face band. Compared with the bulk, both the two sur- 
faces are found to increase the total DOS at the Fermi 
level. The superconductivity of the MgB2 surfaces is 
found to be somewhat strengthened. The work function 
is predicted to be 5.95 and 4.25 eV for the B- and Mg- 
terminated surfaces respectively. The simulated STM 
images are found to be sensitive to the tip-sample dis- 
tances. For the Mg-terminated surface we even observed 
an image reversal for different tip-sample distances. 

We thank Dr. D.M. Chen of the Rowland Institute 
for Science for many helpful discussions. This work 
was partially supported by the National Project for the 
Development of Key Fundamental Sciences in China 



2 



(G1999075305), by the National Natural Science Foun- 
dation of China, by the Foundation of Ministry of Edu- 
cation of China, and by the Foundation of the Chinese 
Academy of Science. The HPCC, NSC, and SC&CG 
Laboratory of USTC are acknowledged for computational 
facilities. 



[1] J. Nagamatsu et al, Naturc(London) 410, 63(2001). 
[2] S.L. Bud'ko et al, Phys. Rev. Lett 86, 1877(2001). 
[3] D.G.Hinks et al, Nature(London) 411, 457(2001). 
[4] J.M. An and W.E. Pickett, Phys. Rev. Lett. 86, 
4366(2001). 

[5] J. Kortus et al, Phys. Rev. Lett. 86, 4656(2001). 
[6] G. Karapetrov et al, Phys. Rev. Lett. 86, 4374(2001). 
[7] G. RubioBollinger et al, Phys. Rev. Lett. 86, 
5582(2001). 

[8] W.N. Kang et al, Science 292, 1521(2001). 
[9] ,J.P. Perdew et al, Phys. Rev. B 46, 6671(1992). 
[10] G. Kressc and J. Hafner, J. Phys.: Condens. Matter. 6, 
8245(1994). 

[11] G. Kresse and J. Hafner, Phys. Rev. B 49, 14251(1994). 
[12] H.J. Monkhorst and J.D. Pack, Phys. Rev. B 13, 
5188(1976). 

[13] A. Lipp and M. Roder, Z. Anorg. Allg. Chem. 344, 
255(1966). 

[14] H. Schmidt et al, Phys. Rev. B 63, 220504R(2001). 
[15] J. Tersoff and D. Hamann, Phys. Rev. Lett. 50, 
1998(1983). 

FIG. 1. Surface band structure of (a) the B-terminated sur- 
face and (b) the Mg-terminated surface. The projected band 
structure of the bulk MgB2 is shown by shadow areas. 

FIG. 2. Charge density contour maps; (a)the x-z plane 
map of the B-terminated surface, (b) the x-z plane map of 
the Mg-terminated surface, and (c) the x-y plane maps at the 
first surface B layers of the B- and Mg-terminated surfaces 
and at the B layer of the bulk respectively. The contours are 
in the unit of electrons/ A"^. 

FIG. 3. Densities of states (DOS) for the surfaces and 
bulk. (al)-(a3) arc the total DOS of the B-terminated sur- 
face, Mg-terminated surface, and bulk. In (a3), we also plot 
the projected DOS for the B and Mg atomic layers in the 
bulk. (bl)-(b3) are the layer resolved projected DOS of the 
first, second and third (B, Mg) layers from the top surface 
of the B-terminated surface respectively. (b4)-(b6) are the 
counterparts for the Mg-terminated surface. The bulk DOS 
was calculated on a 19 x 19 x 19 Monkhorst-Pack k-point grid, 
while the DOS of slab models was obtained on a 19 x 19 x 2 
k-point grid. 



FIG. 4. Simulated STM images of MgB2 (0001) surfaces 
at -|-0.5eV: (a) B-terminated surface with the tip-sample dis- 
tance (d) of 5 A, (b) B-terminated surface with ii=3A, (c) 
Mg-terminated surface with d=5A, and (d) Mg-terminated 
surface with d=3A. 



TABLE I. Relaxation of MgB2 (0001) surfaces (change 
Adi,j of the interlayer distances in percent of the distance 
in the bulk) and work function ($) in eV. 
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Mg-terminated 
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4.25 
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Figure 1 of Li et al. 




Figure 3 of Li et al. 




Figure 4 of Li et al. 



